Contribution of the second Landau level to the exchange energy of the 
three-dimensional electron gas in a high magnetic field 
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We derive a closed analytical expression for the exchange energy of the three-dimensional inter- 
acting electron gas in strong magnetic fields, which goes beyond the quantum limit (L — 0) by 
explicitly including the effect of the second, L — 1, Landau level and arbitrary spin polarization. 
The inclusion of the L — 1 level brings the fields to which the formula applies closer to the laboratory 
range, as compared to previous expressions, valid only for L — and complete spin polarization. 
We identify, and explain, two distinct regimes, separated by a critical density ric. Below Uc, the per- 
particle exchange energy is lowered by the contribution of L = 1, whereas above it is increased. 
As special cases of our general equation we recover various known, more limited, results for higher 
fields, and identify and correct a few inconsistencies in some of these earlier expressions. 

PACS numbers: 71.10.Ca,71.15.Mb,71.70.Di,97.60. Jd 
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I. INTRODUCTION 

The calculation of the exchange energy of an inter- 
acting Fermi gas in high magnetic fields is a fundamental 
problem of many-body physics, with applications in fields 
as diverse as semiconductor physicSfi;^''^ astrophysics, ^'^'^ 
atomic and molecular physicsi*^!^ and density-functional 
theoryiSili^ia 

In a seminal 1971 paper™ Danz and Glasser (here- 
after DG) calculated the exchange energy, ex, of a three- 
dimensional electron gas in high magnetic fields, by 
means of Green's function techniques. A key result are 
analytical expressions for the dependence of on den- 
sity n and magnetic field B, valid if the electrons are 
fully spin polarized and occupy only the spin-down sub- 
level of the lowest Landau level. In an equally impor- 
tant 1974 paper Banerjee, Constantinescu and Rehakii 
(hereafter BCR) also calculated this exchange energy, 
and obtained a result that looks very similar to that 
of DG. The calculations of DG and BCR provided the 
background for a large body of later work on the ex- 
change and correlation energy of the electron gas in 
high magnetic fi.elds^ i^^'^'^i^'^'^^i^^'^^i^^ and are also fre- 
quently quoted as input for the local-density approxima- 
tion to current-density-functional theory in strong mag- 
netic ficlds^^^^S 

A major limitation of these early calculations is their 
restriction to complete spin polarization and to the lowest 
Landau level, which in three dimensions requires either 
magnetic fields that are beyond what is currently achiev- 
able in the laboratory, or restriction to low-density low- 
effective-mass systems. Here we extend the DG many- 
body calculations to the case of arbitrary spin polariza- 
tion, and include the contribution of the second, L — 1, 
Landau level. As a consequence, the range of magnetic 
fields to which the resulting expression applies is ex- 
tended towards weaker fields, as compared to earlier ex- 



pressions. 

Our formula reveals rather complex behaviour of the 
exchange energy, once higher Landau-levels are included: 
As a function of the density n, Cx first drops with a dis- 
continuous derivative at density n^, corresponding to the 
onset of occupation of L = 1, and then passes through 
two regimes, separated by a critical density ric- Below 
Uc, the per-particle exchange energy is lowered (in mod- 
ulus) by occupation of the L = 1 level, whereas above 
ric it is increased (in modulus). The crossover between 
both regimes corresponds to a local minimum of ex{n) 
at Uc- The physics of the drop and of both regimes can 
be understood in terms of the Landau-level structure. 
For currently achievable fields, ric falls into the metallic 
density regime, and thus should be observable. 

Next, we consider various special cases of our general 
expression, in order to make contact with more restricted 
results previously available in the literature, in particular 
those of DG and BCR. Scrutiny of these earlier expres- 
sions reveals a number of small inconsistencies and mis- 
takes, which we correct on the basis of our more general 
expression. 



II. EXCHANGE ENERGY INCLUDING THE 
SECOND LANDAU LEVEL 

The extension of the DG calculation to include higher 
Landau levels had, up to now, not been achieved in closed 
from. DG, and many other workers ^li^ii^iiiii^ii^iilii^ 

go beyond the quantum limit by rewriting the exchange 
integrals via expansion in infinite series, which cannot 
be resummed, or calculate them numerically. Neither 
approach yields analytical expressions that can be used, 
e.g. in the construction of density functionals for current- 
density-functional theory.-Si^ Motivated by the need for 
analytical expressions for lower fields, including L = 1, 
and by the observation of small inconsistencies in avail- 
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able results for L = (see below), we have recalculated 
the exchange energy along the same lines as in the DG 
calculation for L = 0, but keept all contributions from 
the L = 1 level. The result turns out to permit a closed 
analytical expression which, although lengthy, can be ex- 
pressed in terms of the same special functions and phys- 
ical variables as the original DG formula for L = 0. To 
be concise, we here just present the final result. More 
details on its derivation can be found in Appendix [SI 



Three types of terms contribute: one, Cx , arises ex- 
clusively from L ~ 0, another, arises from L = 1, 
and a third, e'^'^^ from inter-level exchange, involving 
contributions from L = and L = 1. As a function of 
the occupation numbers of the spin up and down sub- 
levels of the L = and L = 1 Landau levels, the final 
result for the per- volume exchange energy can be written 



e.(nJ,ni,n^,ni,S) = ^(^)'^[ei°)K^S) + e^l)K,<,S)+eWK,B^^ (1) 
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Here liJc{B) ~ eB/mc is the cyclotron frequency, 
C = 0.57722 is Euler's constant, GH is the Meijer G 
functionSS and Ei is the exponential integral. In prin- 
ciple, this expression hold for arbitrary values of the g- 
factor, as long as L < 1, but we note that if the free- 
electron value g — 2 is employed, an accidental degener- 
acy between the spin-up subband of L = 1 level and the 
spin-down subband of L = 2 occurs. Restriction to L < 1 
is thus only rigorously possible if either g < 2 {& common 
situation in semiconductors) or if the L — 1 level is fully 



polarized, so that its spin-up subband is empty. 

The condition L < 1 implies a restriction on the al- 
lowable values of density and magnetic field. In terms of 
the magnetic length 1{B) — ^JhcjieB) and the density 
parameter rs(n)ao = [3/(4?™)]-'^/^, where ao — fi^ /{me^) 
is the Bohr radius, this restriction is conveniently written 
as 

— < f (1 + V2) (5) 
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for g = 2, and 
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< 
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for g = 0. These conditions are derived in Appendix iBl 

Regarding the spin dependence, we note that spin up 
and spin down contributions do not mix, i.e., their con- 
tribution to Cx can be evaluated separately. Equations 
(H]) to (|4]) permit one to do this for arbitrary values of the 
occupation numbers n1, and thus also for arbitray spin 
polarizations m oc (n^ — n^), where = X^l^'^'^l- 

As illustration of Eq. Fig. [1] displays the exchange 
energy for a combination of densities and magnetic fields 
for which both L — Q and L = 1 levels contribute, and 
compares it with the (erroneous) use of the L — ex- 
pression alone in the same regime. 

Figure [T] reveals interesting behaviour of Cx that ap- 
pears only once L > is allowed for: Upon adding 
more particles (increasing the density), the per-particle 
exchange energy suddenly drops (in modulus) once the 
L — 1 level starts to be occupied at density n^, and then 
passes through two regimes. Initially, it continuous to 
decrease, while for larger densities it increases (in modu- 
lus) up to values larger than those obtained by allocating 
all particles in L = 0. The first regime is entered with a 
discontinuous derivative in the ex{n) curve, indicating a 
zero-temperature phase transition. The second regime is 
entered via a gradual crossover. Both regimes are sepa- 
rated by a critical density ric, at which ex{n) goes through 
a local minimum. 

This intricate behaviour finds its explanation in the 
Landau-level structure: Due to the differences in the spa- 
tial part of the Landau-level wave functions, intra-level 
exchange integrals are larger than inter-level integrals. 
As long as only L = contributes, the exchange energy 
naturally increases with the number of particles. Once 
particles are allocated in the L — 1 level, too, the re- 
duced spatial overlap of their orbitals with those of the 
particles in L = leads to a reduction of the per-particle 
exchange energy. As the number of particles in i = 1 in- 
creases, their intra-level exchange energy also grows and 
overcompensates the initial drop (as seen in the crossing 
of the full and the dashed curve). 

For g = 2, the initial drop is much more pronounced 
than for g < 2, because for this value occupation of the 
spin-down subband of the L — 1 level occurs simultane- 
ously with that of the energetically degenerate spin-up 
subband of the L = level, so that the exchange energy 
of some of the aditional particles is not only reduced by 
small spatial overlap but strictly zero due to zero spin 
overlap with the spin-down particles already occupying 
L = 0. The resulting larger drop may be so large that 
for fields and densities compatible with L < 1, ex(n) 
does not recover the hypothetical value obtained by fill- 
ing only the L = level, as illustrated in the upper set 
of two curves in Fig. [TJ 

This complex behaviour of ex(n) implies that 
magnetic-field Thomas-Fermi-Dirac theory^iii^ii^ 
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FIG. 1: Full curves: Exchange energy per particle, obtained 
from Eq. ([1} by dividing by the density, as a function of rj^ oc 
n^^^ . Dashed curves: (Erroneous) continuation of the L — Q 
expression into the L = 1 regime. Upper set of two curves: 
<7 = 2 and B = 1.4448 x lO^T. Lower set of two curves: q = 
and B = 1.1474 x 10*r. 



current-density-functional theorjt^ii and magnetic-ficld- 
density-functional theory®:'' should have a much richer 
solution space than have ordinary Thomas-Fermi-Dirac 
theory and spin-density-functional theory, where ex{n) 
is strictly monotonous. The density Ud where the 
discontinuous drop occurs, and the critical density Uc 
where Cxin) goes through its minimum, depend on the 
magnetic field. For sufficiently large fields both fall in 
the metallic density range, i.e., both regimes, as well as 
the crossover between them, should be observable in the 
laboratory. 



III. HIGHER FIELDS AND LIMITING CASES 

In the quantum limit, B is so high that only the L ^ 
level is occupied. In this limit, the exchange energy is 
given by keeping only e^"^ in Eq. ([T]). If the free-electron 
value g — 2 is employed, an accidental degeneracy be- 
tween the spin-up subband of L = level and the spin- 
down subband of L = 1 occurs. Restriction to L = 
is thus only rigorously possible if either 17 < 2 or if the 
system is fully polarized, so that the spin-up subband 
of L = level is empty, too. Moreover, the conditions 
guaranteeing L = are stricter than those limiting occu- 
pation to L < 1, given above, and read 
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for generic g G (—2,2). These conditions are derived in 
Appendix [B] As long as these conditions are satisfied, 
Eq. ([1]) with the contribution from just ei*^^ can be used 
for arbitrary occupation of the up and down subbands, 
i.e. for arbitrary spin polarization. 

The further restriction to full spin polarization, i. e. an 
empty spin-up subband, then leads to 
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which is the result obtained by DG.^*' Here the particle 
density n is equal to the occupation ti'^'Zq and we defined, 
following DG, p{n,B) := Sn'^'h^n'^ /{m^to^). As long as 
no higher levels are occupied and 5 = 2, we also have 
p{n, B) — Atp/'hujci where ep is the Fermi energy. (This 
equality breaks down if other spin states or Landau levels 
are involved, ov g ^ 2.) 

If tpj'hLOc ^ 1, then only the bottom of the lowest 
subband is occupied, and the preceding equation reduces 
to 



e^(n, S) : 

In units in which h 
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1, Eq. ITTI) can be written as 
2.11788 , (12) 



(eB)3/2, 

which was first obtained by DG^*' and can be di- 
rectly compared with the corresponding equation (47) 
of BCR,". 

27re2 
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(eB)3/2 



2.32918 



(13) 

featuring a different numerical value inside the brackets. 
This difference was also noticed in Refs. ITsHT^ . 

For ultra-strong magnetic fields (quantified in Fig. ^ 
or for extremely low densities, that value, be it 2.11788 
or 2.32918, becomes negligible relative to the logarithmic 
term, and the DG and BCR expressions become identical, 
both reducing to 
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where ex/Ry and huJc/Ry denote exchange energy and 
cyclotron energy measured in Rydberg [IRy = e'^ /2aQ — 
13.6eV]. In DG this limit appears as their Eq. (1.1), 
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FIG. 2: Applicability of the four expressions discussed here. 
The curves represent the lowest magnetic fields for which the 
indicated expression is valid, for the free-electron value of m 
and g = 2. 



which differs in two ways from Eq. (jl4p : the numerical 
factor inside the logarithm is 0.282 = 2 x 0.141 and the 
prefactor is 27/16, instead of 27/4. 

Our general expression ([1]) contains all these limits as 
special cases, and we can therefore verify which of the 
conflicting expressions is correct. In the high-field limit 
of the quantum limit, we find that the DG Eq. (fT^ is cor- 
rect, while the BCR Eq. ([13]) is not. In the ultra-strong 
field limit, we find that Eq. is correct, while the DG 
Eq. is not. In this latter case, we suspect that DG 
inadvertedly used Hartree units instead of Rydberg units 
{IH — e^/oo = 2Ry is the atomic unit of energy), but 
denoted them as Rydberg, as this would explain both 
the factor of 2 inside and the factor of 1/4 outside the 
logarithm. 

Figure [5] illustrates the magnetic field and density 
regimes for which each of the above equations is valid. 
The curves represent the lowest magnetic fields for which 
the indicated expression is valid, for the free-electron 
value of m and g — 2. For all B{rs) above the full curve, 
L < 1 and Eq. ^ with arbitrary values of tiq, nj, and 
may be applied. For g ~ 2, n[ must be zero to avoid 
degeneracy with the L = 2 Landau level, not included in 
our formula. For 5 < 2, n| is also arbitrary. 

For values of B{rs) above the dashed curve, Eq. 
with arbitrary values of n — jiq may be applied. 



m 

For 



g ~ 2, tiq must be zero to avoid degeneracy with the 
L = 1 Landau level, not included in the DG formula. 
For g < 2, nj is also arbitrary. The relevant expression 
in this case is our Eq. ([T]), but with only the ei"' term 
kept. 

The restrictions L — and L < 1 are precisely 
defined and easily applied. The restriction p <C 1 
which leads from Eq. (fTO|) to Eq. (fTTjl and the condition 
|ln (n/(eS)3/2) | > 2.11788, which leads from Eq. ^ to 
Eq. (|14p . are less precisely defined, and we simply adopt 
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as validity criterium that p < 0.05 and that 2.11788 is less 
than 5% of the logarithmic term. From these criteria, we 
find that for Bfjs) above the dotted curve, the high-field 
limit pT|) becomes valid. The inset shows the values of 
B{rs) above which the ultra- high-field limit p4)) is valid. 

To generate Figs. [T] and [2] we used the free-electron 
value of the electron mass to, requiring very high mag- 
netic fields to satisfy conditions ([5]) to ([5]). Such high 
fields are not without physical relevance. Continuous 
fields ~ 40T and pulsed fields ~ lO^T can be produced in 
the laboratory, white dwarfs have surface fields of order 
~ lO^T, and neutron stars can have fields in excess of 
10^ — lO^r. For effective masses smaller than the free- 
electron mass, the required fields are considerably lower: 
the use of effective masses m* — rescales all B values 
by B ^ so that a reduction of to by a factor of 10 
allows a reduction of S by a factor of 100. Small effective 
masses and low densities thus bring these fields into the 
laboratory range. 



IV. CONCLUSIONS 

Our key result is Eq. ^ , for the exchange energy of the 
three-dimensional Fermi gas in magnetic fields for which 
both the lowest and the second-lowest Landau level con- 
tribute, and the fermions may have arbitrary spin po- 
larization. This equation predicts the existence of two 
physically distinct regimes, one (entered with a discon- 
tinuous derivative) in which occupation of the L = 1 level 
lowers the per-particle exchange energy, and one (entered 
through a gradual crossover) in which it increases it. For 
high, but achievable, fields both the discontinuous drop 
and the gradual crossover between the two regimes occur 
in the metallic density range, so both regimes should be 
experimentally observable. We predict that similar sep- 
aration in two regimes occurs every time a new Landau 
level is included in the calculation. 

The availability of an analytical expression for the 
exchange energy of the electron gas in high magnetic 
fields facilitates the construction of density functional 
for magnetic-field Thomas-Fermi-Dirac theory ,^'^^i^^ 



current-density-functional theory,-'- and magnetic-field 
density-functional theory,—*^ but the intricate structure 
of ex{n, B) implies that their local approximations must 
display a much more complex behaviour than those of or- 
dinary Thomas-Fermi-Dirac and spin-density-functional 
theory. 

This work was supported by FAPESP and CNPq. We 
thank G. Vignale for useful remarks on an earlier version 
of this manuscript. 

APPENDIX A: DETAILS ON THE DERIVATION 
OF EQS. (HI) TO dH) 

From Eq. (2.1) of DG, and following their procedure 
until their Eq. (2.5), we find the per volume exchange 
energy 



Jo (t^ + r'^y^^ r'j 
5^9 (CD sin{r.Dl)L, (I^r^) | (Al) 
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where Q{x) is the step function, Ll(x) are Laguerre poly- 
nomials, and we defined 

Cl^ ^l-(L+^ hcj, - gtioBa, (A2) 




where /iq = eh/2mc is the Bohr magneton. 

Keeping the contributions from the L — and L — 1 
levels and performing the change of variable 
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Using Laplace Transforms to perform the integral,-- we find 
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The integral over the first term, involving only con- 
tributions from L = 0, can be calculated following 
the procedure of Appendix A of the DG paperJ^ The 
other integrals are calculated by integration by parts and 
Laplace transforms. In particular, the integrals involving 
tan~^ [x) can be calculated by using the result of the first 
integral. 

C£ and can be written in terms of the occupation 
number as C7£ = {2TT'^h^l^ /m)nl^ and Dl = 2TiH^nl. 
The step functions turn out to be unnecessary, because 
the term arrising from L — \, and the mixed term con- 
taining contributions from L = Q and L = I, are auto- 
matically zero when = 0. We thus find the exchange 
energy as given Eqs. H]) - (|4|). More details of the calcu- 
lation are available through Ref. [1^ 

An extension of these analytical calculations to higher 
L values seems extraordinarily cumbersome. To include, 
e.g., L = 2 one must keep one more term in the sum of 
Eq. (|Aip . which gives rise to not just one more intra- level 
integral in Eq. (|A5p , but also to many new inter-level in- 
tegrals. In the general case, for arbitrary L, there is no 
guarantee that the resulting integrals can all be solved 
by Laplace transforms, and reduced to some known spe- 
cial function. The three-dimensional case is, in this re- 
gard, more complicated than the two-dimensional onCf^I 
because of the additional dependence of the single- 
particle energies and the resulting Landau-level disper- 
sion. Numerical work on the general case is under way. 



APPENDIX B: DETAILS ON THE DERIVATION 
OF EQS. (HD TO dHD 

The conditions specifying the fields and densities for 
which restriction to L = or L < 1 is valid depend on the 
Landau-level degeneracy. The single-particle dispersion 
in the spin a sublevel of Landau level L is 



2m 



L 



fiLOc + g^ioBa. (Bl) 



In each subband, the Fermi momentum is related to the 
partial density by kpj^ — 2TT'^Pn'^. The degeneracy of 
these levels depends on the value of g. For the free- 
electron value 5 = 2, the spin- up subband of level L and 



the spin-down subband of level L-\-l are degenerate. For 
g = 0, on the other hand, the spin-up and spin-down 
sublevels of Landau level L are degenerate. We first deal 
with g — 2. 

The restriction to L = implies 

= 0,L=l,i) = e(A;, =0,L = 0,T) (B2) 
>ef = e(4o,L = 0,i). (B3) 

From Eq. (|Bip . and using the fact that for L = and 
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g — 2, riQ = n, we immediately find 

1 



27r4?6 ' 



which is equivalent to Eq. ([7]) of the main text. Similarly, 
the restriction to i < 1 implies 

e(fc, =0,L = 2,i) = e(fc, = 0,L = l,t) (B5) 
>ep = e(4i, 1, i) = e(4o, 0, T) = e(4o' 0' i)- (B6) 

The total density is n = tiq + n J + n| . Solving this set of 
equations for n yields 



n < 



V2 



V27r2/3 ' 



(B7) 



which is equivalent to Eq. ([5]) of the main text. 

The corresponding conditions for 5 = follow in the 
same way. The case g = is occasionally used as a 
methodological device in theoretical work, because it al- 
lows to cleanly separate the effects of spin magnetism 
from those of orbital magnetism. Experimentally, 5 = 
occurs, e.g., in systems studied in the context of 5-factor 
engineering of devices. 

When I5I < 2 and 5 ^ 0, no such degeneracy occurs. 
The restriction to L = then implies 
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where n 
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